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TECHNIQUES  FOR  REAL-TIME  OPERATION  OF  FLOOD 
CONTROL  RESERVOIRS  IN  THE  MERRIMACK  RIVER  BASIN 

by 

Bill  S.  Eichert,1  John  C.  Peters,2  and  Arthur  F.  Pabst3 


INTRODUCTION 


This  paper  contains  a  description  of  the  techniques  that  are  under 
development  at  The  Hydrologic  Engineering  Center  for  providing  decision 
criteria  on  a  real-time  basis  for  operating  the  five  flood  control  reser¬ 
voirs  in  the  Merrimack  River  basin.  Techniques  under  development  Include: 

a.  testing  of  alternative  streamflow  forecasting  models 

b.  application  of  computer  program  HEC-5C,  Simulation  of  Flood 
Control  and  Conservation  Systems,  to  develop  decision-criteria  for  system 
operation  on  a  real-time  basis 

c.  use  of  computer  terminals  to  enable  analysis  of  alternative 
forecast  and/or  decision  criteria  in  both  batch  mode  and  interactive 
applications. 

Each  of  these  techniques  will  be  discussed  following  a  brief  description 
of  characteristics  of  the  basin,  reservoir  system  and  automatic  data 
collection  network. 

CHARACTERISTICS  OF  BASIN  AND  RESERVOIR  SYSTEM 

The  Merrimack  River  basin  is  located  in  east  central  New  England  and 
extends  from  the  White  Mountain  area  of  New  Hampshire  southward  Into  the 
northeast  portion  of  Massachusetts.  The  basin  is  134  miles  long  north  to 
south  and  up  to  68  miles  wide  east  to  west,  the  drainage  area  is  about 
5,000  square  miles,  3,800  of  which  are  in  New  Hampshire. 

The  average  annual  precipitation  for  the  Merrimack  River  basin  varies 
from  about  60  inches  in  the  headwaters  to  about  40  inches  in  the  southern 
section,  with  a  basin  average  of  about  45  inches.  Precipitation  is  fairly 
evenly  distributed  throughout  the  year.  During  winter  months  the  precip¬ 
itation  is  mostly  In  the  form  of  snow,  with  amounts  averaging  70  to  100 
Inches  or  more  in  the  north  to  45  to  60  inches  in  the  southern  areas. 
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Floods  can  occur  during  any  season.  The  two  greatest  basin-wide  floods 
occurred  in  March  1936  and  September  1938.  The  1936  event  resulted  from 
two  periods  of  heavy  rainfall  about  a  week  apart  associated  with  significant 
snowmelt.  The  1938  flood  resulted  from  intense  hurricane  rainfall  which 
occurred  after  a  week  of  almost  continual  rain. 

The  reservoir  system  in  the  Merrimack  basin  consists  of  five  reservoirs, 
all  of  which  are  operated  almost  exclusively  for  flood  control.  Drainage 
areas  and  flood  control  storage  capacities  for  the  reservoirs  are  as  follows 


Drainage  Area 

Flood  Control 

Storage 

Reservoir 

sq.  mi. 

ac.  ft. 

in. 

Franklin  Falls 

1,000 

150,600 

2.8 

Blackwater 

128 

46,000 

6.7 

Hopkinton 

382 

70,100\_ 

f.  K 

Everett 

64 

85,500/ 

MacDowel 1 

44 

12,800 

5.4 

A  schematic  diagram  of  the  reservoir  system  is  shown  in  figure  1.  The 
Hopkinton  and  Everett  reservoirs  are  joined  by  a  canal  to  enable  diversion 
from  Hopkinton  to  Everett  during  flood  events.  Franklin  Falls  reservoir 
has  a  relatively  limited  flood  control  capacity  because  it  was  originally 
anticipated  that  another  reservoir  would  be  developed  upstream. 

Population  centers  in  the  basin  are  located,  for  the  most  part,  along 
the  main  stem  of  the  Merrimack  River.  For  computer  simulation  runs 
described  in  this  paper,  the  locations  of  Franklin  Junction,  Concord, 
Manchester  and  Lowell  were  treated  as  damage  centers  (see  Figure  1).  The 
first  three  centers  are  in  New  Hampshire;  Lowell  is  in  northern  Massachu¬ 
setts.  A  major  proportion  of  the  total  dat^.ges  occurs  at  Lowell.  The 
travel  time  from  Franklin  Falls  to  Lowell  is  about  30  hours. 

PRESENT  REGULATION  PROCEDURES 

Because  all  reservoirs  except  Franklin  Falls  contain  about  6  inches  of 
flood  control  storage,  a  large  flood  or  a  series  of  smaller  ones  can  be 
stored  in  four  of  the  reservoirs  without  spilling.  Consequently,  current 
operation  procedures  require  that  the  outlet  works  of  MacDowell,  Blackwater, 
Hopkinton,  and  Everett  reservoirs  be  shut  off  early  in  a  flood  event. 

Because  of  the  limited  storage  capacity  at  Franklin  Falls,  flows  are 
passed  through  this  reservoir  up  to  the  channel  capacity  of  18,000  cfs. 
Inflows  to  Franklin  Falls  are  estimated  from  reservoir  rate-of-rise 
curves;  future  inflows  are  based  on  flows  at  an  upstream  location  called 
Plymouth.  Once  flows  at  Plymouth  have  peaked,  a  release  rate  for 
Franklin  Falls  is  established  based  on  the  estimated  flood  volume. 
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DATA  COLLECTION  SYSTEM 


The  New  England  Division  has  established  a  comprehensive  data  collection 
network  in  the  Merrimack  River  basin  as  well  as  in  several  other  New  Enqland 
basins.  An  automatic  radio  reporting  network  supplies  information  on  rain¬ 
fall  and  river  stage  directly  to  a  computer  in  the  Control  Center  In  the 
Division  office  in  Waltham,  Massachusetts.  Under  computer- programmed  control, 
reporting  stations  can  be  interrogated  singly  or  as  a  group  at  automatically 
selected  time  intervals  ranging  from  six  hour  to  one  hour  periods  based  on 
the  amount  of  river  flow  or  rainfall.  At  present,  river  stage  is  reported  from 
ten  locations  and  precipitation  from  three  locations.  Also  the  New  Enqland 
Division  is  assessing  the  use  of  orbiting  satellites  for  relaying  information 
from  data  collection  platforms.  Four  of  these  platforms  are  currently  in 
use  in  the  Merrimack  basin. 

FORECASTING  TECHNIQUES 

Streamflow  forecasting  may  be  accomplished  using  a  variety  of  techniques. 

A  relatively  simple  technique  would  involve  relating  the  stage  at  a  downstream 
location  to  the  stage  at  some  earlier  time  at  an  upstream  station.  A  rela¬ 
tively  sophisticated  technique  would  involve  modeling  the  precipitation- 
runoff  process  continuously  Involving  all  aspects  of  the  hydrologic  cycle 
deemed  significant.  Accompanying  this  range  of  forecasting  sophistication 
is  a  range  of  required  data.  A  simple  gage  relationship  requires  only 
stage  or  discharge  as  a  function  of  time.  A  precipitation-runoff  model  may 
require  precipitation,  water  equivalent  of  a  snow  pack,  air  temperature, 
dewpoint,  wind  velocity,  insolation,  albedo,  soil  moisture,  frost  depth  as 
well  as  streamflow  discharge.  Forecasting  by  a  simple  model  is  severely 
limited  in  its  capability  to  provide  information  very  far  into  the  future. 

The  sophisticated  precipitation-runoff  model  may  utilize  forecasted  meteo¬ 
rological  conditions  and  provide  forecasts  of  runoff  as  far  as  is  reasonably 
possible  Into  the  future. 

The  selection  of  a  particular  technique  will  depend  on  a  realistic 
assessment  of  the  actual  data  available,  the  accuracies  of  the  forecast 
method,  and  the  use  that  will  be  made  of  the  forecasted  streamflows. 

Streamflow  Extrapolation  Forecasting 

A  method  for  making  relatively  short-term  forecasts  that  utilizes  a 
minimum  amount  of  data  would  be  a  useful  tool  to  complement  a  more  complete 
precipitation  runoff  model.  Initial  efforts  were  directed  to  development  of 
a  simple  forecast  tool  that  would  require  only  observed  streamflow  at  various 
locations  in  the  basin.  The  technique  developed  is  termed  Streamflow  Extrap¬ 
olation  Forecasting.  In  essesnce  it  is  no  more  than  taking  observed  flows  at 
stream  gaging  locations,  extending  them  into  the  future  by  gage  relationships, 
recessing  the  flows  from  that  point  on,  and  routing  them  down  the  basin. 


This  technique  will  be  described  by  use  of  a  simple  illustration. 

Given  only  the  observed  flows  at  stream  gage  locations  A,  B,  C  and  D, 
(Figure  2a)  up  until  the  current  time,  the  problem  is  to  estimate  the 
future  flows  at  sites  A  and  B. 

The  following  series  of  steps  would  be  taken: 

(1)  Extend  the  current  flows  at  A,  3-hours  into  the  future  based  on 
the  immediately  preceding  flows  at  A  and  nearby  station  C.  (Figure  2b) 

(2)  Recess  this  hydrograph  at  A  for  flows  beyond  3  hours  into  the 
future.  This  will  provide  the  forecasted  flows  at  A.  (Figure  2b) 

(3)  Route  the  forecasted  flows  at  A  down  to  B.  (Figure  2c) 

(4)  Extend  the  current  flows  at  B,  6-hours  into  the  future  based 

on  the  immediately  preceding  flows  at  B  and  nearby  station  D.  (Figure  2c) 

(5)  Subtract  the  routed  flows  from  A  from  the  extended  flows  at  B 
yielding  Incremental  local  flows  between  A  and  B.  (Figure  2d) 

(6)  Recess  the  Incremental  local  flow  beyond  6-hours  into  the 
future.  This  will  be  the  forecasted  incremental  flows  between  A  and  B. 
(Figure  2d) 

(7)  Add  these  incremental  flows  to  the  flows  routed  down  from  above 
yielding  the  forecasted  flows  at  B.  (Figure  2e) 

(8)  Continue  downstream  as  required  using  the  preceding  steps. 

The  gage  relationships  used  in  steps  1  and  4,  to  extend  flows  at 

a  given  station,  should  be  developed  from  historical  data.  Multiple 

regression  analysis  may  readily  be  used  to  establish  these  relationships. 

A  future  flow  at  station  A  may  be  correlated  to  current  and  past  flows 
at  stations  A  and  C.  The  time  span  for  extending  flows  (i.e.,  3  hours, 

6  hours)  will  depend  on  the  size,  shape  and  other  characteristics  of  the 
basins. 

Such  a  procedure  which  assumes  little  or  no  future  precipitation  or 

snowmelt  input,  yields  the  minimum  future  flows  in  the  river  system.  If 

future  additional  precipitation  or  snowmelt  occur  the  actual  flows  will 
exceed  the  forecasted  flews. 

Such  a  forecast  can  provide  a  firm  basis  for  establishing  that  a 
reservoir  will  surely  fill.  It  may  not  give  a  long  lead  time  as  to 
when  a  reservoir  will  fill;  but  it  is  free  from  the  uncertainties  of 
having  to  assess  average  basin  precipitation  with  sparse  data,  calculations 
of  snowmelt,  or  of  establishing  losses  for  frozen  or  partly  frozen  ground. 
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The  forecasted  discharges  are  interfaced  to  the  reservoir  operation 
model  through  a  data  file.  This  allows  the  forecast  model  to  run 
Independently  of  the  operation  model.  As  a  separate  program  the  fore¬ 
cast  can  be  made  once  and  then  several  operation  policies  may  be 
evaluated  using  the  discharges.  It  is  not  necessary  to  fit  both  the 
forecast  and  operation  models  in  computer  core  at  the  same  time,  or 
to  resort  to  overlays.  Any  other  forecast  technique  may  be  used  by 
simply  having  the  alternative  model  write  the  appropriately  formatted 
discharge  file. 

Evaluation  of  Forecasts 


In  order  to  choose  between  alternative  forecasting  techniques  it  is 
necessary  to  establish  a  measure  of  forecast  accuracy.  The  measure 
should  reflect  the  error  over  a  certain  finite  period  into  the  future, 
recognizing  that  the  error  over  such  a  span  will  change  as  one  proceeds 
through  the  event.  In  addition  it  would  be  desirable  to  aggregate  the 
error  over  a  series  of  several  floods  rather  than  to  accept  a  given 
technique  only  on  its  reproduction  of  one  event. 

When  historical  data  are  available  the  Streamflow  Extrapolation 
Forecasting  program  will  provide  information  on  the  relative  error 
(eq.  1)  and  the  standard  error  (eq.  2)  over  a  future  time  span 
selected  by  the  user 


i+zl  |Qj  -  Q0BS.| 

j=i+l  J  0OBS.  J 
RELERR1  =  - - - i— 


ST0ERRi 


where 


'll  (Qj  -  QOBSj)2 
3=1+1 


Q  is  forecast  discharge 
QOBS  Is  observed  discharge 
i  is  current  time  index 
i  is  future  span  length  of  forecast 
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Such  errors  can  be  evaluated  at  each  time  period  as  one  steps  through 
the  event.  The  relative  errors  for  the  1936  flood  for  the  Merrimack  River 
at  Lowell,  MA  for  two  forecast  conditions  are  shown  In  figure  3.  Beneath 
the  actual  cumulative  local  flow  hydrograph  are  shown  the  errors  In  a 
forecast  based  only  on  recessed  flows  at  upstream  locations,  and  the 
errors  when  station  flows  are  first  extended  Into  the  future  by  gage 
relationships  and  then  recessed.  Relatively  good  agreement  was  obtained 
for  short  term  forecasts  using  this  technique. 

Another  indication  of  the  adequacy  of  the  forecast  can  be  measured 
in  terms  of  the  efficiency  of  system  operation  by  comparing  the  results 
of  operation  of  historical  floods  with  and  without  forecasts.  One 
measure  of  this  efficiency  that  may  be  useful  for  flood  control  operation 
Is  a  comparison  of  expected  (or  average)  annual  flood  damages  (AAO) 
using  historical  flows  and  forecasted  flows.  Such  a  comparison  Is  made 
In  Table  2  where  the  average  annual  damage  (AAD)  for  run  F-l  Is  about 
the  same  for  locations  8,  9,  and  10  as  the  runs  using  measured  streamflow 
(runs  J  and  17).  The  AAD  for  location  11  is  considerably  higher  for  run 
F-l  than  for  either  run  J  or  17  Indicating  additional  improvement  In 
forecasting  procedures  Is  needed  for  the  location  that  is  substantially 
further  downstream. 
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RESERVOIR  OPERATIONAL  MODEL 
Basic  Objectives  of  Model 

The  HEC-5C  program  was  initially  developed  to  assist  in  planning 
studies  required  for  the  evaluation  of  proposed  changes  to  a  system  and 
to  assist  in  sizing  the  system  components  for  flood  control  and  conser¬ 
vation  requirements.  However,  the  program  can  also  be  used  in  studies 
made  immediately  after  a  flood  to  calculate  the  preproject  conditions 
and  to  show  the  effects  of  existing  and/or  proposed  reservoirs  on  flows 
and  damages  in  the  system.  Special  features  have  been  added  to  the 
program  to  make  it  useful  for  real-time  applications.  The  program  logic 
is  designed  to  minimize  flooding  as  much  as  possible  and  yet  empty  the 
system  as  quickly  as  possible  while  maintaining  the  proper  balance  of 
flood  control  storage  among  the  reservoirs. 

The  above  objectives  are  accomplished  by  simulating  the  sequential 
operation  of  various  system  components  of  any  configuration  for  short 
interval  historical  or  synthetic  floods  or  for  long  duration  nonflood 
periods,  or  for  combinations  of  the  two.  Specifically  the  program  may 
be  used  to  determine: 

a.  Releases  from  reservoirs  during  flood  emergencies  based  on 
local  flow  forecasts  furnished  to  the  program. 

b.  The  evaluation  of  operational  criteria  for  both  flood  control 
and  conservation  for  a  system  of  reservoirs. 

c.  The  influence  of  a  system  of  reservoirs,  or  other  structures  on 
the  spatial  and  temporal  distribution  of  runoff  in  a  basin. 

d.  The  expected  (or  average)  annual  flood  damages  (AAD),  system 
costs,  and  excess  flood  benefits  over  costs. 

e.  Flood  control  and  conservation  (including  hydropower)  storage 
requirements  of  each  reservoir  in  the  system. 

f.  The  determination  of  the  system  of  existing  and  proposed  reservoirs 
or  other  structural  or  nonstructural  alternatives  that  results  in  the 
maximum  net  benefit  for  flood  control  for  the  system  by  making  simulation 
runs  for  selected  alternative  systems. 
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The  input  data  requirements  for  any  basin 'for  HEC-5C  can  be  minimal 
for  preliminary  planning  studies  or  detailed  for  model inq  existing  systems. 
The  minimum  data  requirements0 a re  as  follows: 

a.  General  Information  (4  cards} 

(1)  Title  cards  for  job  (3  cards) 

(2)  Six  miscellaneous  items  including  the  number  of  periods 
of  flow  data,  time  interval  of  flows,  etc. 

b.  Reservoir  Data  (4  cards  per  reservoir) 

(1)  Reservoir  capacities  for  top  of  conservation  and  top  of 
flood  control  elevations. 

(2)  Downstream  control  points  for  which  reservoir  is  to  be 
operated. 

(3)  Reservoir  storage/outflow  tables. 

c.  Control  Point  (including  reservoirs)  Data  (3  cards  per  control 
point) 

(1)  Identification  number  and  title 

(2)  Operating  channel  capacity 

(3)  Channel  routing  criteria  (Muskingum,  modified  Puls,  Working 
R/D,  Tatum,  or  Straddle-Stagger) 

d.  Flow  Data 

Inflow  or  local  flow  data  (or  observed  flows  and  reservoir  releases) 
for  each  control  point  for  one  or  more  historical  (including  forecasted 
flows)  or  synthetic  floods. 

Optional  data  on  flood  damages  may  also  be  used  by  inputting  peak 
discharge-damage  data  where  flood  damages  are  directly  related  to  maximum 
stage  (or  discharge)  obtained  during  a  flood  event. 

General  Operational  Criteria  for  Model 

a.  Reservoirs  are  operated  to  satisfy  constraints  at  individual 
reservoirs,  to  maintain  specified  flows  at  downstream  control  points,  and 
to  keep  the  system  in  balance.  Constraints  at  individual  reservoirs  are 
as  follows: 
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(1)  When  the  level  of  a  reservoir  is  between  the  top  of  conser¬ 
vation  pool  and  the  top  of  flood  pool,  releases  are  made  to  attempt  to 
draw  the  reservoir  to  the  top  of  conservation  pool  without  exceeding  the 
designated  channel  capacity  at  the  reservoir  or  at  downstream  control 
points  for  which  the  reservoir  is  being  operated. 

(2)  Releases  are  made  equal  to  or  greater  than  the  minimum 
desired  flows  when  the  reservoir  storage  is  greater  than  the  top  of  buffer 
storage,  and  or  equal  to  the  required  flow  if  between  level  one  (top  of 
inactive  pool)  and  the  top  of  buffer  pool.  No  releases  are  made  when  the 
reservoir  is  below  level  one.  Releases  calculated  for  hydropower  require¬ 
ments*  will  override  minimum  flows  if  they  are  greater  than  the  controlling 
desired  or  required  flows. 

(3)  Releases  are  made  equal  to  or  less  than  the  designated  channel 
capacity  at  the  reservoir  until  the  top  of  flood  pool  is  exceeded,  then 

all  excess  flood  water  is  dumped  if  sufficient  outlet  capacity  is  available. 

If  insufficient  capacity  exists,  a  surcharge  routing  is  made.  Input  options 
permit  channel  capacity  releases  (or  greater)  to  be  made  prior  to  the  time 
that  the  reservoir  level  reaches  the  top  of  the  flood  pool  if  forecasted 
inflows  are  excessive. 

(4)  The  reservoir  release  is  never  greater  (or  less)  than  the 
previous  period  release  plus  (or  minus)  a  percentage  of  the  channel 
capacity  at  the  dam  site  unless  the  reservoir  is  in  surcharge  operation. 

b.  Operational  criteria  for  specified  downstream  control  points  are 
as  follows: 

(1)  Releases  are  not  generally  made  (as  long  as  flood  storage  remains) 
which  would  contribute  to  flooding  at  one  or  more  specified  downstream 
locations  during  a  predetermined  number  of  future  periods  except  to 

satisfy  minimum  flow  and  rate-of-change  of  release  criteria.  The  number 
of  future  periods  considered  is  the  lesser  of  the  number  of  reservoir 
release  routing  coefficients  or  the  number  of  local  flow  forecast  periods 
specified  on  input  data. 

(2)  Releases  are  made,  where  possible,  to  exactly  maintain  down¬ 
stream  flows  at  channel  capacity  (for  flood  operation)  or  for  minimum 
desired  or  required  flows  (for  conservation  operation).  In  making  a 
release  determination,  local  (intervening  area)  flows  can  be  multiplied 
by  a  contingency  allowance  (greater  than  1  for  flood  control  and  less 
than  1  for  conservation)  to  account  for  uncertainty  in  forecasting  these 
flows. 


*No  Corps  hydropower  projects  are  in  the  Merrimack  River  Basin. 
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c.  Operational  criteria  for  keeping  a  reservoir  system  ir>  balance 
are  as  follows: 

(1)  Where  two  or  more  reservoirs  are  in  parallel  operation  above 
a  common  control  point,  the  reservoir  that  is  at  the  highest  index  level, 
assuming  no  releases  for  the  current  time  period,  will  be  operated  first 
to  try  to  increase  the  flows  in  the  downstream  channel  to  the  target  flow. 
Then  the  remaining  reservoirs  will  be  operated  in  a  priority  established 
by  index  levels  to  attempt  to  fill  any  remaining  space  in  the  downstream 
channel  without  causing  flooding  during  any  of  a  specified  number  of  future 
periods. 


(2)  If  one  of  two  parallel  reservoirs  has  one  or  more  reservoirs 
upstream  whose  storage  should  be  considered  in  determining  the  priority 
of  releases  from  the  two  parallel  reservoirs,  then  an  equivalent  index 
level  is  determined  for  the  tandem  reservoirs  based  on  the  combined 
storage  in  the  tandem  reservoirs. 

(3)  If  two  reservoirs  are  in  tandem  (one  above  the  other),  the 
upstream  reservoir  can  be  operated  for  control  points  between  the  two 
reservoirs.  In  addition,  when  the  downstream  reservoir  is  being  operated 
for  control  points,  an  attempt  is  made  to  bring  the  upper  reservoir  to 
the  same  index  level  as  the  lower  reservoir  based  on  index  levels  at  the 
end  of  the  previous  time  period. 

Use  of  Contingency  Allowance  and  Foresight 

Two  key  input  items  are  used  in  determining  reservoir  releases  based 
on  downstream  flooding  as  discussed  under  "operational  criteria  for  down¬ 
stream  control  points."  These  factors  are  the  number  of  future  time 
periods  (IFCAST)  that  should  be  checked  for  possible  future  flooding 
(called  forecast  periods)  and  the  contingency  allowance  (CFLOD)  which 
is  multiplied  times  the  cumulative  uncontrolled  downstream  flow  to 
account  for  uncertainty  in  forecasts.  For  simulation  of  historical 
floods  (where  flows  are  known  for  duration  of  flood)  a  contingency 
factor  of  1  and  an  infinite  forecast  period  could  be  used  in  order  to 
operate  with  maximum  foresight.  However,  these  assumptions  would  not 
simulate  "real  world"  conditions  where  large  errors  in  forecasting 
future  streamflows  are  possible.  These  two  key  factors,  for  the 
simulation  of  historical  floods,  should  be  selected  so  that  the  oper¬ 
ational  efficiencies  in  the  planning  mode  will  approach  the  expected 
efficiencies  under  flood  emergency  conditions.  During  flood  emergencies 
these  factors  should  be  used  to  insure  that  the  forecasting  errors  do 
not  cause  reservoir  releases  to  be  made  which  will  cause  major  unneces¬ 
sary  flood  damages.  The  sensitivity  of  the  system  to  different  values 
of  these  two  factors  can  be  determined  by  simulating  the  operation  and 
resulting  flood  damages  for  a  series  of  different  sized  flood  events  for 
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the  system.  The  difference  between  the  average  annual  damages  (AAD)  for 
various  combinations  of  these  factors  will  help  to  evaluate  the  sensitivity 
of  these  factors.  Table  1  illustrates  how  the  reservoir  system  responds  to 
these  factors.  The  adopted  value  for  the  number  of  forecast  periods  (IFCAST) 
was  four  and  the  adopted  contingency  factor  (CFLOD)  was  1.2.  Because  a  time 
interval  of  3  hours  was  used,  the  duration  of  the  adopted  forecast  period  was 
12  hours. 


TABLE  1 

AVERAGE  ANNUAL  DAMAGE  VS  FORECAST  PERIOD 
AND  CONTINGENCY  FACTOR  USING  HISTORICAL  FLOWS 

AAD 


RUN 

CFLOD 

IFCAST 

(in  $1000) 

1 

1.0 

2 

1191 

2 

1.0 

4 

1083 

3 

1.0 

6 

1048 

4 

1.0 

10 

1070* 

5 

1.2 

2 

1116 

6 

1.2 

4 

1046 

7 

1.2 

6 

1032 

8 

1.2 

10 

1052* 

9 

1.2 

20 

1099* 

10 

1.4 

2 

1082 

11 

1.4 

4 

1047 

12 

1.4 

6 

1056* 

13 

1.4 

10 

1090* 

14 

1.6 

4 

1075 

♦One  would  expect  these  values  to  be  less  than 
the  previous  values.  They  are  not  because  for 
the  larger  events  a  long  forecast  period  causes 
reservoir  releases  to  be  diminished  relatively 
early  in  the  event.  When  the  reservoirs  event¬ 
ually  go  uncontrolled,  the  resulting  flooding 
is  greater  than  would  have  occurred  if  the 
releases  had  not  been  diminished  early  in  the 
event.  The  increase  in  damages  in  the  larger 
events  exceeds  the  decrease  in  damages  for  the 
smaller  events. 
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TABLE  2  -  AVERAGE  ANNUAL  DAMAGES  SUMMARY  -  FORECASTED  FLOW 


12 


Use  of  Forecasted  Flows 


If  flow  forecast  models  are  available,  t'.e  same  type  of  simulation  runs 
can  be  made  to  determine  the  proper  forecast  period  and  contingency  factor 
by  using  forecasted  streamflow  for  one  or  more  historical  floods  and  one  or 
more  ratios  of  those  floods  to  calculate  the  average  annual  damages.  In 
most  cases,  the  best  operation  should  occur  where  the  average  annual  damages 
are  a  minimum.  Table  2  illustrates  results  using  forecasted  flows.  The 
adopted  values  for  the  forecast  period  (IFCAST)  were  four  3-hour  periods 
and  contingency  factors  (CFLOD)  were  assumed  as  1.2  and  1.4  respectively. 
Results  indicate  forecast  flows  are  generally  adequate  for  locations  8,  9, 
and  10,  but  not  at  location  11. 

The  adopted  values  for  historical  floods  (where  future  flows  are 
known)  should  not  necessarily  be  the  same  as  the  adopted  values  during 

flood  emergencies  since  the  forecasted  flows  will  not  be  the  same  as  the 

observed  historical  flows.  In  general,  the  contingency  factor  for  historical 
floods  should  be  selected  to  produce  the  same  AAD  as  the  run  using  fore¬ 
casted  flows.  The  number  of  periods  of  foresight  should  be  the  same 

regardless  of  the  source  of  flows. 

OPERATIONAL  CRITERIA  FOR  MERRIMACK  BASIN 

An  essential  task  associated  with  computer  simulation  of  the  Merrimack 
reservoir  system  is  evaluation  of  input  parameters  for  computer  program 
HEC-5C  to  obtain  the  most  desirable  operation  of  the  system.  Some  of  the 
key  input  parameters  are  listed  in  Table  3.  Alternative  operation  criteria 
were  evaluated  by  determining  average  annual  damages  based  on  spatial  and 
temporal  runoff  variations  associated  with  the  March  1936  and  September 
1938  flood  events.  While  average  annual  damage  is  a  useful  criteria  for 
selecting  operating  policies,  other  factors  such  as  legal  and  political 
considerations  must  also  be  used  in  the  evaluation.  The  procedure 
used  in  HEC-5C  for  estimation  of  average  annual  damages  is  as  follows: 

a.  Ratios  are  determined  for  application  to  selected  historical  flood 
events  that  are  representative  of  the  full  range  of  frequency  of  flood 
occurrence.  For  example,  ratios  of  1.4,  1.0,  0.8,  0.7  and  0.5  were  applied 
to  inflow  and  local  flow  hydrographs  for  the  March  1936  flood  to  obtain 
five  floods  for  which  system  operation  was  to  be  simulated.  Frequencies 
associated  with  peak  discharges  for  the  five  floods  were  determined  from 
frequency  curves  for  unregulated  flows  for  locations  where  damages  were 

to  be  computed. 

b.  Reservoir  system  operation  is  simulated  for  each  flood,  that  is, 
for  each  set  of  inflow  and  local  flow  hydrographs  obtained  by  applying 
ratios  to  hydrographs  for  a  historical  event.  Frequencies  associated 
with  peak  "regulated"  discharges  are  assumed  to  be  the  same  as  the 
"unregulated"  frequencies.  Figure  4  illustrates  natural  and  regulated 
frequency  curves  for  Lowell.  Points  on  the  regulated  frequency  curve 

in  figure  4  represent  peak  discharges  resulting  from  system  simulation 
for  a  specific  set  of  operation  criteria. 
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c.  Damage-discharge  relations  input  to  the  computer  program  enable 
determination  of  dollar  damages  corresponding  to  the  peak  discharges  at 
each  damage  center  for  each  flood.  Figure  5  illustrates  the  damage- 
discharge  relation  that  was  used  for  Lowell.  This  is  an  approximate 
relation  that  will  be  updated  in  the  future. 

d.  Damage-frequency  relations  are  established  for  each  damage  center 
for  both  natural  and  regulated  conditions.  Figure  6  illustrates  these 
relations  for  Lowell.  The  computer  program  integrates  the  area  below 
the  damage- frequency  curves  to  obtain  average  annual  damages. 

The  average  annual  damage  calculation  is  influenced  by  the  distribution 
of  runoff  for  the  historical  "pattern"  storms.  Some  characteristics  of 
runoff  production  for  the  1936  and  1938  floods  can  be  ascertained  from  the 
hydrographs  in  Figures  7  and  8.  These  plots  show  discharge  per  square  mile 
for  inflow  to  Franklin  Falls  reservoir,  unregulated  flow  on  the  Contoocook 
river  at  Penacook  and  uncontrolled  local  flow  (runoff  from  all  areas  down¬ 
stream  from  nearest  upstream  reservoirs)  at  Lowell. 

The  1936  flood  reflects  high  runoff  production  over  the  entire 
Merrimack  basin.  The  flatness  of  the  peak  for  local  flow  at  Lowell 
reflects  the  relatively  slow  responsiveness  of  this  portion  of  the  basin. 
Lowell  is  a  key  location  because  a  large  proportion  of  total  damages 
occurs  there.  Consequently,  an  objective  in  operating  the  reservoir 
system  is  to  try  to  avoid  "building  on"  the  local  peak  at  Lowell. 

Figure  8  indicates  that  runoff  production  from  the  Contoocook  was 
relatively  high  for  the  1938  event.  However,  this  portion  of  the  Merrimack 
basin  is  'controlled'  with  four  of  the  five  reservoirs. 
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TABLE  3 


INPUT  PARAMETERS  FOR  HEC-5C 

1.  Number  of  periods  of  future  (forecasted)  flows  that  will  be 
used  to  determine  reservoir  releases. 

2.  Contingency  factors 

These  are  ratios  to  be  applied  to  flows  in  determining 
reservoir  releases;  factors  are  used  to  account  for  limited 
knowledge  of  future  flows  beyond  the  forecast  period. 

3.  Control  points  for  which  reservoirs  are  to  be  operated. 

4.  Rate-of-change-of-release  criteria  for  reservoirs. 

5.  Channel  capacity  criteria  for  control  points. 

6.  Minimum  release  vs  reservoir  elevation  criteria. 

7.  Pre-release  criteria 

a.  Whether  or  not  pre-releases  will  be  permitted.  A  pre¬ 
release  is  a  flood-producing  reservoir  release  that  is  made 
when  the  reservoir  level  is  below  the  top  of  flood  control 
pool.  The  release  is  based  on  the  anticipated  flood  volume 
exceeding  available  capacity. 

b.  Reservoir  elevation  that  pre-releases  will  be  geared  to. 

c.  Magnitude  of  pre-release  permitted  (can  be  specified 
as  a  function  of  reservoir  elevation). 
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RESULTS  OF  AVERAGE  ANNUAL  DAMAGE  RUNS 

In  setting  up  average  annual  damage  runs  for  HEC-5C,  a  variety  of 
approaches  could  be  used  in  selecting  floods  and  flood  ratios.  For  example, 
one  or  more  ratios  of  a  number  of  different  historical  events  could  be 
incorporated  in  a  single  average  annual  damage  computer  run.  Another 
approach  is  to  determine  average  annual  damages  for  separate  sets  of  ratios 
applied  to  individual  historical  events.  Results  of  these  individual 
average  annual  damage  runs  could  be  weighted  depending  on  how  representative 
individual  storms  are  of  the  overall  flood-producing  characteristics  of  the 
basin.  Two  separate  sets  of  simulation  runs  were  made  to  compute  average 
annual  damages  for  the  Merrimack  basin.  As  indicated  previously,  one  set 
is  based  on  using  five  ratios  of  the  1936  event.  A  second  set  uses  five 
ratios  of  the  1938  event. 
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Results  of  the  average  annual  damage  runs  are  summarized  in  Table  4. 

The  "base"  runs,  labeled  A  in  Tahle  4,  were  made  with  HEC-5C  input 
parameters  specified  as  follows. 

a.  A  forecast  period  of  12  hours;  that  is,  discharges  up  to  12  hours 
in  the  future  were  considered  in  determining  reservoir  releases. 

b.  A  contingency  factor  of  1.2  was  applied  to  local  flows  for  purposes 
of  reservoir  release  determination. 

c.  The  reservoir  system  was  operated  for  all  control  points  shown  in 
figure  1 . 

d.  Rate-of-change  of  release  criteria  were  specified  so  as  not  to  be 
a  constraint  on  releases  from  Franklin  Falls  reservoir. 

e.  Fixed  channel  capacities  were  specified  for  all  control  points 
based  on  information  supplied  by  the  New  England  Division. 

f.  A  table  of  values  for  minimum  permissible  release  as  a  function 

of  reservoir  elevation  was  specified  for  Franklin  Falls  reservoir  as  shown 
on  Table  4  (note  2) . 

g.  Pre-releases  were  permitted;  at  Franklin  Falls  reservoir,  pre¬ 
releases  were  made  if  inflows  to  the  reservoir  during  the  12-hour  forecast 
period  would  cause  the  reservoir  level  to  rise  above  elevation  394  (5  feet 
above  the  spillway  crest). 

The  discharge-damage  relation  for  Lowell  (Figure  5)  that  was  input  to 
HEC-5C  had  a  maximum  discharge  ordinate  of  180,000  cfs.  Because  discharges 
larger  than  180,000  cfs  were  used  in  the  damage  analysis,  the  discharge- 
damage  relation  was  extrapolated  by  the  computer  program  as  shown  in 
Figure  5.  The  effect  of  using  the  alternative  extrapolation,  also  shown 
in  Figure  5,  on  average  annual  damages  was  less  than  21  for  both  natural 
and  regulated  flows. 


18 


Paper  12 


Operation  criteria  for  average  annual  damage  runs  other  than  the  base 
runs  are  summarized  in  Table  5.  Some  observations  pertaining  to  results  of 
the  average  annual  damage  simulation  runs  are  as  follows: 

a.  Average  annual  damages  based  on  floods  patterned  after  the  March  1936 
flood  are  of  approximately  the  same  magnitude  as  average  annual  damages  based 
on  floods  patterned  after  the  September  1938  flood. 

b.  Operational  criteria  used  for  the  base  run  (run  A)  produced  the 
lowest  average  annual  damages  for  floods  patterned  after  the  March  1936 
flood;  operational  criteria  that  does  not  utilize  the  pre-release  option 
(run  E)  produced  the  lowest  average  annual  damages  for  floods  patterned 
after  the  September  1938  flood  and  only  slightly  more  damage  for  the  1936 
f 1 ood . 


c.  Of  the  order  of  75ft  of  the  total  average  annual  damages  occurs  at 
Lowell  on  the  basis  of  the  approximate  stage-damage  relationship  for  that 
location. 

d.  Damages  associated  with  very  large  floods  account  for  a  major 
proportion  of  average  annual  damages;  this  is  illustrated  in  figure  9  which 
shows  the  relation  between  percent  of  average  annual  damage  and  recurrence 
interval  at  Lowell  for  the  base  run  for  floods  patterned  after  the  March 
1936  flood,  (e.g.,  45ft  of  average  annual  damages  occur  under  regulated 
conditions  from  floods  having  a  recurrence  interval  of  300  years  or  greater). 

e.  Significant  discharge  and  damage  reduction  at  the  1000  year  flood 
level  is  due  to  the  surcharge  storage  available  in  the  reservoirs  due  to 
the  limited  discharge  capacity  of  the  uncontrolled  spillways  since  the 
flood  control  storages  were  exceeded  very  early  in  the  largest  floods. 
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TABLE  5 


OPERATION  CRITERIA  FOR  SIMULATION  RUNS 

Cri teria 


Same  as  for  base  run,  except  minimum  release  not 
specified  for  Franklin  Falls  reservoir. 

Same  as  for  base  run,  except  relatively  large 
minimum  releases  were  specified  (as  a  function  of 
reservoir  elevation)  for  Franklin  Falls  reservoir 
(see  note  2  of  Table  4  for  values). 

Same  as  for  base  run,  except  relatively  small 
minimum  releases  were  specified  (as  a  function  of 
reservoir  elevation)  for  Franklin  Falls  reservoir. 

Same  as  base  run,  except  pre-releases  were  not  made. 

Same  as  base  run,  except  reservoir  system  was  not 
operated  for  Lowel 1 . 

Same  as  base  run,  except  pre-releases  were  not 
made  and  minimum  releases  for  Franklin  Falls 
reservoir  were  not  specified. 

Same  as  base  run,  except  pre-releases  were  not 
made  and  reservoir  system  was  not  operated  for 
Lowell . 

Same  as  base  run,  except  reservoir  elevation  of 
389  was  used  at  Franklin  Falls  reservoir  for  pre¬ 
release  determination. 

Same  as  base  run,  except  minimum  releases  of  18,000 
cfs  from  Franklin  Falls  reservoir  were  specified, 
system  was  not  operated  for  Lowell. 

Same  as  base  run,  except  the  pre-release  option  was 
modified  to  include  volume  of  recession  of  hydrograph 
past  the  period  of  forecast. 
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OUTPUT  DISPLAYS 


r 


Batch  Mode  vs  Interactive  Mode 


Execution  of  a  computer  program  in  batch  mode  requires  that  all  input 
for  a  computer  run  be  supplied  to  the  computer  prior  to  program  execution. 
An  interactive-mode  execution  is  where  the  user  can  interact  with  the 
computer  during  the  execution  of  a  job.  As  used  in  the  application 
described  herein,  the  interactive  mode  is  used  to  selectively  print  out 
data  from  an  output  file  of  the  system  operation  that  has  been  generated 
in  batch  mode.  This  enables  the  user  to  review  any  portion  of  the  output 
that  he  desires.  The  output  file  can  be  permanently  saved  and  interrogated 
at  future  times  from  one  or  more  computer  terminal  sites. 

While  output  displays  from  high-speed  line  printers  used  in  the  batch 
mode  can  provide  any  amount  of  output  desired,  the  level  of  output  must 
be  specified  prior  to  making  the  computer  run.  Presently,  after  a  run  has 
been  made,  output  not  previously  requested  can  only  be  obtained  by  making 
another  complete  simulation  run.  An  alternative  method  would  be  to  save 
the  output  file  and  print  in  batch  mode  by  writing  a  special  program.  If 
the  turn-around  time  is  adequate  (say  less  than  30  minutes)  and  the  program 
execution  cost  is  small,  then  the  batch  mode  is  the  best  way  to  get  the 
necessary  output  assuming  that  a  high-speed  printer  is  available.  Where 
batch  mode  turn-around  times  are  long,  or  high-speed  printers  are  not 
readily  available,  the  slow-speed  terminals  can  be  an  effective  way  of 
obtaining  a  limited  amount  of  information  rapidly.  After  looking  at 
selected  data  through  the  slow-speed  terminal  the  output  can  then  be 
directed  to  a  high-speed  printer  if  desired. 


I 
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The  subroutine  PROUT  in  HEC-bC  is  used  to  print  output  for  the  high¬ 
speed  printer  (batch  mode),  slow-speed  teletype  terminal  (interactive  mooe) 
and  cathode  ray  tube  terminals  (interactive  mode).  All  output  devices  can 
be  used  to  print  any  combination  of  types  of  output  (see  samples  on 
Figures  10-11)  except  for  the  graphical  plots  (see  Figure  17)  available 
with  cathode  ray  tube  terminals.  Printer  plots  (Figure  12)  can  be  requested 
by  the  batch  mode  printers.  The  types  of  output  that  can  be  requested  are 
as  follows: 

OUTPUT  DESCRIPTION 

*  Input  Card  Listing 

*  Input  Flows 

*  Input  Data  for  System  Specification 

*  Output  -  Normal  Sequential  (by  control  point) 

*  Output  -  Reservoirs  -  by  Period 

*  Output  -  Reservoir  Releases  -  by  Period 

*  Output  -  Reservoir  Regulation  Summary  -  Single  Flood 

*  Output  -  Reservoir  Regulation  Summary  -  All  Floods 

*  Output  -  Hydrologic  Efficiencies 

*  Output  -  Computer  Check  for  Possible  Errors 

*  Batch  Economic  Summary 

*  User  Designed  Output  -  Results  by  Period 

*  User  Designed  Output  -  Summary 

Interactive  Terminal 

While  the  batch  mode  selects  desired  output  by  input  cards,  the  slow- 
speed  terminal  asks  the  user  questions,  as  illustrated  in  figure  13, 
concerning  what  type  of  function  (see  Figure  14)  should  be  performed 
next  (output  from  operation,  modify  input  deck  HEC-5C,  etc.).  If  oper¬ 
ational  output  is  desired,  the  type  of  data  (option  number  of  Figure  14) 
and  possibly  variable  codes  and  locations,  on  Figure  15,  as  well  as  the 
output  mode  (plot,  tabulate  or  save  on  tape)  are  required  to  be  specified 
to  the  computer.  Data  can  be  input  to  the  computer  through  the  terminal 
by  depressing  the  appropriate  terminal  keys  or  (for  certain  CRT's)  by 
touching  an  electronic  pen  to  the  appropriate  instruction  (see  Figure  15) 
on  a  selection  list  (menu)  which  lies  on  a  graphic  tablet.  After  the 
tabulated  output  (see  Figure  16)  or  plotted  results  (Figure  17)  are 
obtained,  additional  data  can  also  be  selected  and  printed  as  desired. 

Any  combination  of  data  types  (reservoir  outflow,  storage,  elevation, 
downstream  flow,  etc.)  for  any  location  can  be  tabulated  or  plotted. 
Depending  on  the  size  of  the  paper  (or  screen)  up  to  10  different  items 
can  be  tabulated  side  by  side  and  up  to  5  different  time  dependent 
variables  can  be  plotted  on  a  single  graph.  Two  different  scales  can 
be  used  for  the  plots  as  shown  on  Figure  17  where  the  inflow  and  outflow 
are  plotted  on  the  left  scale  (discharge)  and  the  reservoir  level  on  the 
right  scale.  When  all  of  the  desired  output  has  been  displayed,  the  data 
can  be  transferred  to  the  line  printer  at  a  nearby  batch  location. 
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In  the  batch  mode,  data  cards  for  toe  forecast  program  are  read  by  the 
card  reader  and  an  output  file  of  the  forecasted  basin -wide  flows  is 
obtained.  The  HEC-5C  data  cards  can  be  loaded  at  the  same  time  or  at  a 
later  time,  and  the  system  simulation  is  performed  for  the  duration  of 
the  forecasts  and  the  output  is  directed  to  the  line  printer.  Any 
desired  change  in  the  forecast  or  operation  requires  a  few  new  data 
cards  and  rerunning  either  the  operation  model  or  both  the  forecast 
and  operation  models.  Where  adequate  computer  turn-around  is  available 
this  process  can  be  accomplished  in  30  minutes  or  less. 

With  a  slow-speed  terminal  the  same  cycle  can  be  accomplished  by 
using  the  keyboard  or  the  electronic  pen.  In  addition,  high  quality 
graphical  displays  of  the  status  of  the  system  can  be  obtained.  After 
the  operational  data  has  been  displayed,  function  5  of  figure  14  can  be 
selected  and  an  interactive  program  called  REVISE  will  allow  the  revision 
and/or  execution  of  the  data  files  of  either  the  forecast  or  operation 
models.  It  is  then  possible  to  once  again  display  selected  output.  A 
diagram  of  this  process  is  shown  as  figure  18. 

COMPUTER  SYSTEMS  REQUIRED  FOR  REAL-TIME  OPERATION 

In  order  to  reap  the  benefits  of  the  real-time  operation  tools  described 
previously,  access  to  digital  computer  equipment  is  a  necessity.  This 
access  may  take  several  different  forms  depending  on  characteristics  of 
the  data  acquisition  system,  the  forecasting  and  operation  models  used, 
and  available  communications  equipment.  Real-time  water  resource  operations 
may  make  use  of  only  in-house  equipment,  only  remote  site  equipment,  or 
some  combination  of  each. 

The  first  basic  function  to  be  accomplished  is  that  of  acquisition  of 
available  data.  Information  required  may  include  observer4  collected  data 
reported  by  voice  communications  and  analog  or  digital  signals  received  by 
appropriate  equipment.  Such  information  may  be  received  over  dial-up 
phone  lines,  dedicated  phone  lines,  radio  repeating  links,  or  satellite 
repeating  links. 

When  the  volume  of  data  is  great  it  is  obviously  desireable  to  have 
the  data  recorded  directly  on  a  medium  that  may  be  read  by  machines.  Even 
more  desireable  is  to  have  the  whole  data  acquisition  function  occur  under 
control  of  a  mini-computer.  If  this  is  done  the  mini -computer  may  also 
perform  several  functions  such  as,  error  checking,  data  reduction,  permanent 
logging,  report  generation,  etc.  In  most  cases  a  mini-computer  dedicated 
to  such  functions  would  be  required. 
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Once  the  data  in  reduced  form  (i.e.,  discharge,  precipitation,  etc.) 
is  available,  the  forecast  may  be  performed.  Computer  equipment  for  this 
function  will  vary  depending  on  the  size  of  the  forecast  program.  Some 
forecasting  models  may  execute  well  on  the  same  mini -computer  used  for 
data  acquisition.  Others  would  require  such  extensive  reprogramming  in 
order  to  operate  on  a  mini -computer  that  other  alternatives  would  be 
desirable.  In  some  cases  additional  large  scale  in-house  computers  may 
be  available.  In  most  cases  such  facilities  will  not  be  directly 
available  to  district  offices.  The  use  of  larger  capacity  remote  site 
comput4ng  becomes  very  attractive  in  such  a  situation.  A  portion  of  the 
reduced  data  may  be  sent  to  the  remote  site  directly  by  the  mini-computer, 
or  read  in  from  paper  tape,  or  magnetic  cassette.  The  forecast  may  then 
be  performed  at  the  remote  site  with  results  returned  to  the  local  site 
and/or  saved  for  future  reference  at  the  remote  site. 

A  similar  situation  exists  for  executing  large  system  operation  programs. 
Again  remote  site  computing  offers  a  cost  effective  solution.  The  forecasted 
flows  may  be  passed  on  to  the  operation  model  through  common  access  to  a  data 
file.  Results  from  the  operation  routine  may  then  be  returned  to  the  local 
site,  or  be  held  for  future  reference.  Display  of  the  results  may  be 
performed  by  the  in-house  mini-computer,  or  be  handled  directly  from  the 
remote  site. 

Problems  of  using  a  remote  site  computing  center  around  two  factors; 

(1)  guaranteed  access  to  facilities  twenty-four  hours  per  day,  three  hundred 
and  sixty-five  days  per  year,  and  (2)  reliable  communications  and  power 
supply  even  under  extremely  adverse  conditions.  It  is  impossible  to 
guarantee  access  to  any  one  facility  at  all  times.  It  is  possible  however, 
to  have  access  to  several  remote  facilities  and  thus  provide  as  high  a 
level  of  "guaranteed"  access  as  deemed  necessary.  Communications  and 
power  supply  which  involve  ground  lines  are  quite  vulnerable  to  inter¬ 
ruption  during  major  storm  activity.  Backup  power  supply  may  be  easily 
supplied  by  emergency  generators. 

Backup  communications  based  on  a  ground  network  may  require  an  individual 
to  physically  carry  necessary  data  to  an  alternate  input  site  outside  the 
area  affected  by  the  communications  interruption.  An  attractive  backup  for 
communications  which  is  rapidly  being  developed  is  to  communicate  to  the 
remote  site  by  way  of  a  satellite  link. 

The  use  of  large  scale  machines  for  the  forecasting  and  operations 
aspects  of  real-time  operation  are  attractive  because  of  the  ease  of  up¬ 
dating  or  improving  the  models  used.  Other  modeling  techniques  may  be 
quickly  compared  and  substituted  for  those  currently  in  use.  When  mini¬ 
computers  are  used  for  executing  large  programs,  changes  to  the  program 
often  entail  major  restructuring  of  overlays. 
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FUTURE  WORK 


It  is  planned  to  implement  the  forecast-operation-display  capabilities 
described  in  this  paper  at  the  Control  Center  of  the  New  England  Division 
in  the  near  future.  The  next  step  will  be  to  interface  these  capabilities 
with  the  existing  automated  data  collection  system. 

Alternative  forecasting  techniques  other  than  the  streamflow  extra¬ 
polation  procedure  described  herein  will  be  tested.  Application  of  the 
computer  program  Streamflow  Synthesis  and  Reservoir  Regulation  (SSARR) 
developed  by  the  North  Pacific  Division  is  anticipated. 

The  entire  procedure  for  real-time  simulation  will  be  thoroughly 
tested  and  "fine-tuned"  once  it  is  operational  at  the  Control  Center. 
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FIGURE  8 
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INTERACTIVE  HEC-5C  MENU 

if  you  want  data  by  period  FUNCTION  2  (USER  DESIGNED  OUTPUT) 
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INTERACTIVE  HEC-5C  OUTPUT 
SAMPLE  -  USER  DESIGNED  TABULATION 
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